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The universal conductance fluctuations of quasi-two-dimensional systems are analyzed with ex- 
perimental considerations in mind. The traditional statistical metrics of these fluctuations (such as 
variance) are shown to have large statistical errors in such systems. An alternative characteristic 
is identified, the inflection point of the correlation function in magnetic field, which is shown to 
be significantly more useful as an experimental metric and to give a more robust measure of phase 
coherence. 
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Universal conductance fluctuations (UCF) are the 
static but sample-dependent variations in electrical con- 
ductance that arise in diffusive mesoscopic systems due 
to quantum interference EHU Certain aspects of UCF are 
highly sensitive to phase coherence, and so the analysis 
of UCF is an important experimental tool in determining 
the decoherence resulting from a dynamic environment, 
which is a key factor in any quantum device. Although 
there are other transport signatures of coherence, they ei- 
ther require a specific sample geometry (e.g., Aharonov- 
Bohm rings 4 ) or fail to provide coherence information in 
environments lacking time-reversal symmetry (this is the 
case with weak localization, for example^). The statis- 
tics of UCF, on the other hand, provide a general-purpose 
measure of coherence for any mesoscopic diffusive device. 

In the past, experimental studies of UCF have primar- 
ily been carried out with quasi-one-dimensional (quasi- 
1D) samples, i.e. wires or channels that are very narrow 
compared to the coherence length. The theory of UCF is 
well developed for this caseP^I Unfortunately, the quasi- 
1D regime has several practical drawbacks that have lim- 
ited its applicability as a tool to study coherence. First, 
the typical scale of conductance fluctuations in magnetic 
field is large (hundreds of millitesla in some cases) and 
hence very large field ranges are required for a proper sta- 
tistical analysisP In semiconducting systems, the useful 
field range is often limited by the onset of Landau quanti- 
zation, so it may be impossible to gather sufficient statis- 
tics for a reliable measurement P The requirement of large 
fields also obscures features of the coherence time that 
may be field dependent.^ Another issue is that the sys- 
tem must be very narrow to avoid a complicated crossover 
to the quasi-two-dimensional (quasi-2D) regime, so edge 
effects can become important. 

The quasi-2D regime, in which the coherence length 
is smaller than the device lateral dimensions, does not 
share the disadvantages of the quasi- ID regime men- 
tioned above. Although the fundamental theory of UCF 
has been well-established for over 20 yearsf^ to this day 
its consequences for the quasi-2D case are only partially 
resolvedP 

The goal of this manuscript is to present a detailed 
analysis of the quasi-2D UCF correlation function, with 
an eye toward using this statistical analysis in experi- 



ments. It is shown that phase coherence information can 
be reliably extracted only from the correlation function 
with respect to magnetic field. The traditional metrics 
of this correlation function — variance and half-width — 
are shown to have unexpectedly large statistical errors, 
which limit their utility in an experiment. The inflection 
point of the magnetic field correlation function, on the 
other hand, turns out to be a particularly robust mea- 
sure of coherence, which extends the practical use of UCF 
as a coherence detector to quasi-2D systems. The appli- 
cability of UCF as a thermometeiQ23 in the quasi-2D case 
is also discussed. 

Section [I] identifies several experimentally-accessible 
metrics of the quasi-2D UCF correlation function and 
describes their dependence on decoherence and temper- 
ature. Section [n] investigates the statistical errors that 
naturally arise in the analysis of UCF data, and how 
those errors affect the accuracy of various techniques for 
measuring coherence from UCF correlations. Section |III| 
discusses the adjustment of the correlation function that 
occurs when there are static symmetry-breaking interac- 
tions, such as spin-orbit interaction. Section [TV| compares 
these results to the quasi- ID case. 

I. CORRELATION LENGTHS OF QUASI-2D 
CONDUCTANCE FLUCTUATIONS 

A typical UCF experiment begins with measuring con- 
ductance G(/x, -B), which is a function of the externally 
controlled chemical potential /i and the magnetic field 
B perpendicular to the sample plane. The ensemble- 
averaged conductance background G is subtracted off to 
yield the conductance fluctuations SG(fi, B) = G(/z, B) — 
G(/j,,B). The fluctuations contain a large amount of in- 
formation, but most of it is sample specific and chaoti- 
cally sensitive to the exact disorder configuration. The 
statistics of 6G, on the other hand, are not sample spe- 
cific and are encapsulated by its correlation function^ 

F(8n, SB) = IgJJjl, B) 8G{^ + 6/j,,B + SB). (1) 

which is independent of \x and B under usual conditions, 
but does depend on temperature T, dephasing rate r^" 1 , 
and diffusion constant DP^l 
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TABLE I. Asymptotic field and energy correlation lengths. 
Prefactors have been numerically determ ined , whereas expo- 
nents are analytically derived (see Appx. A 3 1. 



FIG. 1. Dephasing rate dependence of several characteristic 
scales of the quasi-2D UCF correlation function in magnetic 
field. Results for inflection point SBi (solid blue), roundness 
8B r (dotted red), and half- width SB i (dashed black) are in- 
dicated by thick lines. Thin lines indicate the asymptotic 
forms in Table U) The inset shows the correlation function for 
k B T — Shr^ 1 and graphically depicts the definitions of SBi , 
5B r , and SBi. 



The theoretical two-parameter correlation function can 
be computed using various methods (see Refs. 11121 and 
Appx. [A]) , and in principle can be compared directly to 
its experimental counterpart in order to extract tem- 
perature or phase coherence information. In practice, 
however, experimentalists usually analyze one-parameter 
cross sections for convenience: F^Sfi) = F(Sfi,0), or 
Fb{SB) = F(0, SB). We first focus on the magnetic field 
cross section, F B {8B). 

It is customary in UCF studies to characterize Fb(SB) 
by its half- width, and compare to values provided by the- 
ory. In the quasi-2D case, however, different field scales 
associated with the correlation function depend in differ- 
ent ways on t^ 1 , and it is therefore important to identify 
which scale is most appropriate for a particular experi- 
ment. We consider three different field scales of the cor- 
relation function (Fig. 1 inset): 

• The half-width SB i, defined by F B (SBi) = 

^F B (0), is the point where correlation has fallen 
to 50% of the variance. 

• The roundness 5B r = \ 2F B (0)/F£ (0) | s charac- 
terizes correlations at very small field separation, 
where F'^(SB) = d 2 F B /dSB 2 . 

• The inflection point SBi, defined as the point where 
F^(SBi) — 0, is the field separation at which cor- 
relation falls the fastest (the minimum F' B (8B)). 
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Figure [T] shows how these three field scales, calculated 
from the theoretical F B (5B), depend on t^ 1 in the case 
that only one dephasing rate is relevant (the case of mul- 
tiple dephasing rates is discussed in section III). Imme- 



diately one can see that the three scales are not propor- 
tional, illustrating the multi-scale nature of Fb(SB). The 
field scales are expressed here in terms of a characteris- 
tic thermal field k B T/ (2eD), and the dephasing rate in 
terms of the thermal time, to make this plot general for 
any quasi-2D system. 

Table [I] lists the asymptotic behaviour of each field 
scale in the thermally smeared limit (hr^ 1 <C k B T) and 
the unsmeared limit {hr^ 1 ^> k B T). The unsmeared 
limit is rarely encountered at low temperatures, where 
dephasing is typically dominated by the contribution of 
electron-electron interactions^ giving t^ 1 — aksT/h, 
for some a less than unity. In the smeared limit, SBi 
depends equally on T and (cf . Ref . [5| , while round- 
ness SB r depends logarithmically on T. Remarkably, SBi 
has no direct T-dependence in either limit. This is a de- 
sirable characteristic because a measurement of SBi then 
yields the value of r^ 1 directly, without needing exact 
knowledge of T. 

Energy correlation lengths {5/J.i , Sfi r , S/ii} are shown 
in Fig. [2j computed from the theoretical F /1 (Sfi), follow- 
ing definitions analogous to the SB correlations. Asymp- 
totic forms are listed in Table [TJ For strong thermal 
smearing (k B T ^> frr^ 1 ), _F /i (<5/i)/i ; ) t (0) approaches a 
universal function of S/i/ '(fcgT), which is independent of 
t7 ; this gives the universal correlation lengths S/j,i ri 
listed in Table [I] As a result, F /1 (S^) is not useful for mea- 
suring t^ 1 , but can instead be used as a thermometer.^ 

The strategy of using UCF as a primary thermometer 
has been shown to be effective for quasi-lD systems.^ 
For the quasi-2D correlation function, however, conver- 
gence to the universal form is very gradual, and all three 
metrics deviate significantly from their asymptotic values 
even when hr^ 1 < 0.1k B T (Fig. 2k. The deviation is par- 
ticularly severe for the metric ojti identified in Ref. I10[ 

e.g., 20% at hr^ 1 = 0.05fc B T. Somewhat faster conver- 
gence is observed for the inflection point 6 fa (e.g., 8% 
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FIG. 2. Mapping of several energy correlation lengths to 
dephasing rate, analogous to Fig. [I] Results for 8/ii (solid 
blue), 8/i r (dotted red), and 6fxi (dashed black) are indi- 
cated by thick lines. Thin lines indicate the asymptotic 
forms in Table [I] The inset shows the correlation function 
for ksT = Shr^ 1 and graphically depicts the definitions of 
Sfii, 8fi r , and Sfii. 



deviation at Ht^ 1 = O.OSfcsT). For the highest accuracy, 
both 8Bi and should be measured: together these 
provide unique values for both T and t&. 



II. STATISTICAL ERRORS IN QUASI-2D UCF 
MEASUREMENTS 



Statistical errors play a major role in experimental 
studies of UCF. Even when the conductance G(fi,B) is 
measured exactly (without noise), there are two statisti- 
cal error sources that affect the quality of the analysis: 
random errors due to a limited data set, and system- 
atic errors due to the background subtraction procedure. 
A formal treatment of these errors can be found in Ap- 
pendix [b] 

In this section we examine the two types of errors as 
they apply to experiments measuring Fb(SB) from quasi- 
2D UCF in magnetic field. It is worth noting that er- 
rors associated with measurements of F^d^i) are very 
different, but this function contains little phase coher- 
ence information. A careful statistical analysis is re- 
quired to estimate the errors in Fb(SB) with any de- 
gree of accuracy — a nonintuitive result that comes from 
the multi-range nature of the function. Inspection alone 
generally overestimates the number of effectively inde- 
pendent samples that contribute to an averaged correla- 
tion function, often by orders of magnitude. Moreover, 
this number depends strongly on which field scale is be- 
ing extracted from the correlation function (5Bi vs. SBi, 
etc.). 
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FIG. 3. Example of statistical errors in UCF analysis for the 
case r=lK,T^ = 100 ps, and D = 0.03 m 2 /s. The solid line 
indicates the UCF correlation function Fb(SB). The dotted 
lines show the autocorrelations of two simulated G(B) traces 
over a 1 T range in B, with mean background subtraction 
(Btot = -Bsm = 1 T; see Appx. 
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The dashed line is an av- 
erage of 100 such autocorrelations (B t ot = 100 T, B sul — 1 T). 
Inset: Expansion of the boxed region, graphically depicting 
the magnitudes of the two types of error in variance. 
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FIG. 4. Guidelines for the required total scan length in field 
(Btot) for reaching 10% standard deviation in the various mea- 
sures of Fb(SB). As expected, the standard deviation of any 



measure falls as B t , 



when Btot is increased. 



Random errors appear in the correlation function when 
it is estimated from a finite data set. Figure [3] shows an 
example of random errors, which appear as fluctuations 
in the estimated correlation function. These fluctuations 
in turn cause uncertainties in all derived parameters such 
as correlations lengths or variance. The magnitudes of 
the random errors depend on the total scanned range of 
data, i?tot, which may be distributed over multiple inde- 
pendent scans. As seen in Fig. |4j the total scan length 
required for a reliable estimate is different by many or- 
ders of magnitude for the various statistical metrics. 
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FIG. 5. Guidelines for the required smoothing-length in field 
(-B am ) for reaching —10% bias in the various measures of 
Fb(SB). Here we have taken mean background subtraction 
with a scan length B aln . Biases in Bi, B r , Fb(0) fall as 
w log(_B sm )/_B sm when B BU1 is increased. Bias in Bi falls as 
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FIG. 6. Contributions to the magnetic field correlation 
function from three UCF modes at T = 2 K, D = 
0.03 m 2 /s. The dotted curves show separate modes with 
t~I s = {20 ns" 1 , 220 ns" 1 , 2120 ns" 1 } with iVi, 2 , 3 = {1,1,2}, 
and the solid curve is their sum [Eq. |2|]. 



III. INFLUENCE OF SYMMETRY-BREAKING 
DISORDER 



Systematic errors (biases) occur when the background 
conductance is estimated from measured conductance 
data itself by fitting, smoothing, etc. Background sub- 
traction inevitably affects the correlation function: any 
smoothed background removes slowly-varying but gen- 
uine conductance fluctuations that are longer than a 
"smoothing length" B sm , which must be determined from 
the smoothing algorithm and is necessarily less than B to t 
(see Appx. B 1 1 . Considering the effects of a smoothed 



background, the consistent loss of long-ranged fluctua- 
tions distorts the analysis in the form of a downward bias 
(see, e.g., Fig. [3]). This bias directly impacts the accu- 
racy of not only variance (Fb(0)) but also SBi and SB r , 
whose definitions rely on variance (Fig. [5]). SBi is nearly 
immune to biases since it does not depend on variance. 



Figures 3-5 demonstrate the dramatically different 
sensitivity to errors for half- width and inflection point. 
To put this difference in practical terms, consider a typ- 
ical low temperature UCF measurement of in a dis- 
ordered semiconductor, where one might have T = 1 K, 
T = 100 ps, and D = 0.03 m 2 /s. To achieve 10% ac- 
curacy (systematic error) using SBi , an extremely large 
B S ui — 3 T would be required even though the value 
of SBi itself is just 5.2 mT. Using SBi, on the other 

hand, 10% accuracy would be obtained for B sm = 5 mT, 
smaller by three orders of magnitude compared to the 
SBi case. Similarly, 10% precision in using SBi would 
require a total scan length B to t — 300 mT, compared to 
B tot = 45 T for SBi. 



Many systems of practical interest involve more than 
one dephasing rate, giving a UCF correlation function 
that deviates from the behaviour described in Sec. |U 
These additional dephasing rates arise from symmetry- 
breaking static disorders that partially dephase UCF, 
e.g., frozen magnetic i mpuri ties spin-orbit coupling, 7 
or valley-orbit coupling P3EU With such disorder, the cor- 
relation function becomes the sum of independent modes 
with a set of distinct dephasing rates {t" 1 } and degen- 
eracies {Ni}, in the form 



F(S/j,,5B) 



mFlr^KS^SB) 
+ N 2 F[r^ 1 ](S f i,SB) 
+ N 3 F[r^}(8fi,SB) 



(2) 



These modes are oft en k nown as the diffuson singlets 
and diffuson triplets EEMH1 With this summation, the 
field scales of F are determined by a complicated mix- 
ture of the temperature and various dephasing rates, so 
that the considerations of Sec. [I] may not directly apply. 
Nevertheless, each independent mode F[ri] is of the type 
described in Appendix|Aj so it is straightforward to com- 
pute Eq. pi) and numerically extract the field scales. 

Figure [6]shows an example involving the three modes 
appropriate for graphene. We have chosen typical val- 
ues for T — 2 K: a decoherence rate of t^ 1 = 20 ns" 1 , 
an intervalley rate of t iv = 100 ns" 1 , and an intravalley 
rate of r* = 2000 ns _1 . The resulting modal dephas- 
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and 



mg rates are 

t^ 1 = t^ 1 + + t" 1 . Although the rates are greatly 
different in magnitude, each mode has a significant con- 
tribution to Fb(SB) because of thermal smearing (see 
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FIG. 7. Effect on field correlation's inflection point from the 
combination of two UCF modes, in smeared (fcsT > Whr^ 1 ) 
and unsmeared (ksT < 0.02/iTj -1 ) limits. Dashed lines show 
the unperturbed inflection point when there is no secondary 
mode (N 2 = 0). 



Appx. A3). This causes the variance, half- width, and 



roundness to differ greatly from the value expected of the 
dominant mode (ti) alone. For instance, SBi is twice as 
large, which would be misinterpreted (by the considera- 
tions of Sec.|T]) as a dephasing rate four times larger than 
the actual rj~ . The inflection point, however, remains 
a reliable measure of r-f 1 even when the additional rates 
are neglected, with only a 4% error. 

Figure [7] examines how SBi , computed for the case of 
only two modes, depends on the relative dephasing rates 
of the two modes. Here t^ 1 might be the decoherence 
rate from dynamic scatterers that affects all UCF modes, 
while (t^ 1 — Tj - ) could be the extra static symmetry- 
breaking rate affecting only the second mode. When the 
symmetry-breaking rate is comparable to the decoher- 
ence rate, the inflection point may be displaced from the 
value expected of r-f alone. The degree of displacement, 
however, never exceeds a factor of 1.8 (this only occurs if 
N 2 /Ni = 3 and hr^ 1 <C fc^T). Fig. [7] also gives a simple 
rule for SB^. the secondary mode can be neglected when 
it dephases much more rapidly than the primary mode, 
t^ 1 > lOrf 1 ; such a simple rule does not apply for other 
aspects, e.g., Fb{0) or SBi. In the opposite limit, when 

dephasing rates for all modes are similar, t^ 1 ~ T-j -1 , the 
field scales are determined by the considerations of Sec. [I] 
with the single dephasing rate r" 1 (or r ^" 1 ), and Fg(0) 
is only increased by a trivial factorP^EH 

The statistical errors (Sec. [TTJ) are also modified by the 
presence of symmetry-breaking disorder, often becoming 
much larger. The considerations of Appx.|B]may be used 
to evaluate errors in a general correlation function with 
multiple dephasing rates. 



IV. COMPARISON TO QUASI-1D CASE 

The characteristics of quasi-2D UCF can be compared 
to the quasi- ID regime, which occurs when the mate- 
rial is shaped as a long and very narrow strip with a 
width W that is much smaller than the dephasing length 
Ltf, = \J -Dt^, yet where the length between contacts is 
longer than L$. An advantage of quasi- ID UCF (espe- 
cially useful in metals) is that shaping the material into a 
wire produces a lower background conductance which al- 
lows UCF to appear with higher contrast. The quasi- ID 
correlation function Fg(5B), which is known exactly in 
both the smeared and unsmeared limits, demonstrates es- 
sentially single-scale behaviour, falling as 1/|(5B| 3 at high 
SB^This explains why the half- width performs well as a 
measure of coherence in the quasi- ID system: its statis- 
tical errors are only somewhat higher than the inflection 
point for a given range of field. 

For completeness, we note the values of inflection point 
for the quasi-lD system in the 'dirty' regime, where the 
elastic mean free path is smaller than W. The quasi- ID 
energy correlation function in Ref. [10] gives the inflec- 
tion point 5[ii = O.549/i/t0 in the unsmeared limit; this 
converges on the universal value Sfii — 2.14/csT in the 
smeared limit. As for the the magnetic field inflection 
point, the formulas in Ref. [6] yield SBi = \fihj(eWL$) 
in the unsmeared limit, and SBi = \eW L^) in the 
smeared limit. 



V. CONCLUSION 

In closing, the analysis of UCF is a valuable tool for 
experimental studies of dephasing, as it can provide in- 
formation on both the electron temperature and the de- 
phasing time under a wide range of conditions. We have 
shown that the correlation function inflection points in 
energy and magnetic field provide an accurate and effi- 
cient way to do so. This inflection point method will be 
especially useful in the experimental study of the mag- 
netic field dependence of coherence in quasi-2D systems. 



Appendix A: Theoretical quasi-2D UCF correlation 

The procedure for calculating the theoretical F was 
originally formulated by Stone, Lee, and Fukuyama^ 
and first computed for the quasi-2D case by Bergmann. 2 
We present a different numerical approach than that in 
Ref. [21 by first calculating the unsmeared correlation 
function analytically. Of special interest is the asymp- 
totic behaviour in the smeared limit, discussed in sub- 
section IA 31 

The main result is ( A4 ) , which when combined with 



the thermal smearing convolution ( A7 ) provides the cor- 
relation function of UCF (A8|. The correlation function 



here is computed for the case of a single dephasing rate, 
and may be used as a building block in more complicated 



G 



situations of multiple dephasing rates (Sec. Ill) or other 
modifications (e.g., Zeeman splitting). 



1. Exact unsmeared correlation function 

Consider a rectangular conductor with length L x be- 
tween the source and drain contacts (at opposing edges) 
and width L y . We begin by analyzing the fluctuations 
in the quantity Gq(E,B), which is the source-drain con- 
ductance at a fixed energy E (with thermal smearing not 
yet taken into account) and in magnetic field Bz. The 
ensemble-averaged correlations in the conductance fluc- 
tuations, SGq = Go — Go, are given by^ 



F (8E, SB) = SG (E, B) SG (E 
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e 4 4£> 2 v 
h 2 Li ^ 
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where G is a constant prcfactor depending on the number 
of intact symmetries in the system.^ The A„ are the 
eigenvalues of the diffusion equation 



D ( iV - -A(r) 
h 



5E 



Qif) = A„Q(rO, 



(A2) 

where Vxi = SBz, and with Dirichlet and Neumann 
boundary conditions at the contact and vacuum edges, 
respectively. 

The primary assumption of the quasi-2D case is that 
we can ignore the boundary conditions, in which case the 
eigenfunctions Q(r) of (A2| are Landau levels! 3 ® The 
"cyclotron frequency" here is 2De\SB\/h, giving the se- 
ries of eigenvalues 



A, 



2De\SB\ 



SE 



for non-negative integers n, with the degeneracy of 
e\5B\L x L y /h for each level. This places L x and L y de- 
pendence only in the prcfactor of Fq, allowing us to write 

F (SE,SB) = C^ L ^Ue,(3). 

where we have introduced a scale-independent function 
$o(e,P) in terms of a dimensionless energy and field, 

e = SE- t^/H, = \SB\ ■ 2eD T4> /h. 



The sum in ( Al I may be solved exactly in terms of 



the complex digamma function 'ijj(z) and its derivative, 
ip'(z). For So this solution is written as 
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This expression does not evaluate when either P or e are 
zero; taking limits, one obtains the variance 3o(0, 0) = ^- 
and one-parameter correlational! 
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2. Thermal smearing effect 

The direct effect of temperature is to average the bare 
conductance Go over a range of values SE ~ S/j, ± fcsT. 
In detail, the measured conductance G is determined by 
the thermal smearing convolution 

/oo 
dEf F (E-l*)G (E,B), 
-00 

where f' F (SE) = 4fc ^ r sech 2 ( 2fc ^ T (5.E) is the Fermi func- 
tion. The smeared correlation function, F(Sfi,SB), 
is then obtained by a convolution of F (SE, SB) with 
the function -j^r^-^p) where k(x) = |(|coth| — 
l)/sinh 2 f. 3 

In terms of scale-independent variables, the thermal 
smearing effect modifies 5o — ► Si by the convolution: 

$<z{e',p)= [ O °d S K( ^ r) -d (e-e',P). (A7) 
J —00 ^ 

where we have defined a dimensionless temperature 
X = k B T ■ T<j,/h. 



(A3) Equation ( A7 ) is quickly computed by way of fast fourier 



transforms and the convolution theorem!^ The scaled 
result of thermal smearing is then written as 



F(5n,8B) = C 



e 4 LyDrtf, 



d%(e',P). 



(A8) 



h? LI 

Note that here, e' = SfXT^/h (not SE). 

3. Highly smeared behaviour 

We now examine the behaviour of quasi-2D UCF un- 
der a large amount of thermal smearing, I > 1, As 
seen in Eq. |A5|, $ (e,0) ps |2£| _1 for la rge e, and this 
plays a critical role when the convolution ( |A7[ ) is applied. 
The 1/1 e| behaviour is responsible for the long- ranged 
and multi-scale nature of Fb(SB) (Sec. [i]), its unusually 
large statistical errors (Sec. Ill]) , and the high sensitivity 
to broken- symmetry modes (Sec. III). 

An immediate consequence of the l/|e| behaviour 
is the logarithmic form of the variance under thermal 
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smearing! 2 ^ The asymptotic behaviour of variance (from 
Eq. (A7|) in the smeared limit is 



ffi(0,0)« — ln(C X), 



T» 1 



(A9) 



where we find numerically that Cq — 4.1. This gives an 
extremely weak r^-dependence of the measured variance, 
F(0,0) ocln(Tr )/T. 



More generally, the correlation function (A4) falls as 
■$o( e i/3) ~ min{|e| — , f3~ 1 } when either argument is 
large. In this case we can approximate (A4) by taking 
| + (1 + ie)/P » | + ie/ /3, giving a r^-independent form: 



5o&P)™20K™/0), /3»lore»l 

where Euler's reflection formula gives f(x) = ^ tanhx + 
\ sech 2 x. When a high amount of thermal smearing is 
applied, the first term in f(x) dominates. This gives the 
intermediate-field behaviour of the field correlation func- 
tion in the smeared limit, 



calculating the inaccuracy of the ergodicity assumption 
under weak violations. 

A typical experiment measures conductance G{x) over 
a limited range x = • • • L. Here a: is a parameter such as 
fx or B. Next, a background estimate gB(x) is computed 
from G(x), then subtracted to yield the estimated fluc- 
tuations, Sg(x) = G(x) — gg (x). Finally, the correlation 
function is estimated as 



f(Sx) 



1 



L — 5x 



L — 5x 



Ax Sg(x)5g(x + Sx), 



and the correlation length estimators xi, x r and Xi are 
extracted from f(Sx). This estimator f(Sx) differs from 
the true correlator F(Sx) — 5G(x) SG(x + Sx) for two 
statistical reasons. First, Sg(x) has lost some information 
from the true fluctuations SG(x) due to the background 
subtraction, which leads to systematic error in f(Sx). 
Second, the limited range L leads to random errors in 
f(Sx) depending on the particular realization of 5G(x). 
These two error mechanisms will be addressed separately. 



62; 



(A10) 



where we find numerically that C\ — 29.1. From (A9) 
and ( A10 ) one obtains the asymptotic half- width listed 
in Table^ /3i » Ci(Z/C )*. 

The correlation function derivative -&^o{s,P) ap- 
proaches zero rapidly, as e sech 2 (7re//3), for large s. As a 
result the thermal smearing convolution imposes a simple 
behaviour for this field derivative: 

±fa(e,/3)*W)?&p-, T»/3andT»l (All) 

for a function [)(/?) = J^^ds A$o(s, /3). This separa- 
ble j3 dependence explains why the inflection point /3j is 
constant for large T. The zero-field curvature can also be 
computed, ^^(0, P)\p=o « t/(0)/(6T), which provides 
the smeared-limit asymptotic roundness when combined 
withEq. (|A9j, f3 r « [21n(C T)/() / (0)] a5 . Numerically we 
find E)'(O) = 1/12. 



Appendix B: Theory of statistical errors in 
measurements of correlation functions 



This appendix explores the statistical errors that oc- 
cur in the estimation of the correlation function from a 
generalized fluctuating quantity G(x). We will use upper- 
case G and F to represent the error-free quantity and its 
ideal correlation, and lower-case g and / to represent es- 
timated values. The overline notation (SG SG, /, x~l etc.) 
in this section refers specifically to an ensemble average 
(average over disorder configurations). The ergodicity 
assumption, which equates this with an average over x, 



1. Background subtraction errors 

To calculate the background subtraction bias, we con- 
sider the simplest possible procedure which is to sub- 
tract the mean of G(x) over the measured interval L. 
The results of this section can be extrapolated to a gen- 
eral background fitting procedure by taking an effective 
L « L sm . For example, if a polynomial fit to G(x) is 
subtracted, then L sm = Ltotai/^, where -Ltotai is the mea- 
sured range and n is the number of degrees of freedom in 
gB(x) [e.g., n — 3 for a parabolic gs(x)]. 

The mean background is 5s(x) = J Q dx' G(x'), 

which gives an error [Sg(x) — 5G(x)] — ~j; Jo & x ' 5G{x > ) 
in the estimated fluctuations. The resulting systematic 
error in the autocorrelation function f(Sx) is approxi- 
mately constant for Sx <C L: To first order, 

XL 



Bias{/(fe)} = f{Sx) - F(Sx) a -f (0)-^. (Bl) 

Here, ij, is a characteristic correlation length, defined 
basically by the area under the correlation function: 



XL 



zhF(z) 



dz 1 
_ L V L I F(0) 



(B2) 



For short-ranged correlation functions, Xl would be a 



constant for large L, and so the systematic error (Bl) 



does not apply in subsection B 1 where we are essentially 



would fall as l/L; this bias then would be similar to the 
well-known sample variance bias from independent sam- 
ple statistics, agreeing with the intuition of G(x) con- 
taining "many independent fluctuations" , each having 
length xl- In the quasi-2D UCF case, however, F(Sx) 
only falls as 1/Sx [see (A5), (A6)] so the value of xl 
diverges logarithmically as L increases. Hence the sim- 
ilarity with independent sample statistics does not hold 
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for the quasi-2D UCF variance bias, as there is no well- 
defined "independence length" . 

The bias in variance /(0) leads to direct effects on 
the half- width estimator (xi) and the roundness esti- 

v 2 ' 

mator (x r ), as these are both sensitive to the abso- 
lute variance. The roundness estimator x r is biased by 
Biasjx,.} « —^x t xl/L. The bias in xi is given by 



Biasja; i } 



F(0) x L 
2F'(xi) L ■ 



The inflection point estimator Xi depends only on the 
second derivative of f(8x), so to first order x% has no bias; 
taking into account higher order terms omitted from (Bl ) 
we obtain Bias{x 4 } w -(2/L 2 )F(x i )/F"'(x i ). 



2. Random errors 

Next we suppose the background has been determined 
perfectly, giving us the exact fluctuations: Sg(x) = 
8G(x). Although we have f(Sx) — F(Sx) in this case, 
the measured f(Sx) will have random deviations from 
F(Sx) due to the limited data set. The random fluctua- 
tions in f(Sx) can be expressed in terms of a two-point 
correlator f{8x\)f(8x2)- If 8G{x) is gaussian (as is the 
case for UCF^) then Isserlis' theorem yields 



f(Sx 1 )f(8x 2 ) - }{8x x ) f(8x 2 ) 

= ~[H(8x 1 -Sx 2 ) + H(8x l + 8x 2 )], (B3) 

for a dataset of large length L 3> 8x±, 8x 2 , where H(88x) 
is a higher-order correlator, defined as 



H{8Sx) 



dzF(z)F(z + SSx). 



(B4) 



Equation ( B3 ) and its derivatives allow the determina- 



tion of random errors in any aspect of f(Sx), including 
its correlation lengths. For instance, the random error in 
variance [/(0)] is given by Var{/(0)} = 2H(0)/L. 



The half- width estimator xi is sensitive to the errors 

2 

in both f(xT) and /(0), modulated by the local slope 
F'(xT), giving 

Var{xi} = Var| ^ ( _ } j 



1 IH(0)+H(2xT)-2H(xT] 



L 



F'{xlf 



(B5) 



The error in the inflection point estimator cc, depends 
on the fluctuation of f"(xl), modulated by the local slope 
F"'{xl). 



Varfo} = Vax{f'(xi)/F'"(xi)} 
_ 1 H""{Q) + H""(2x-) 
~ L F'"{xl) 2 ' 
The roundness estimator, x r — -y/2/(0)/|/"(0)|, is in 
fluenced by changes in both /(0) and /"(0): 



(B6) 



Var{x r } = Var ( - ^ 

IV2F(0)|F"(0)| #W 



H(0) 



F(0)H""(0) 2H"(0) 



F(0)\F»(0)\ 



\F h 



F't (0)2 

' (B7) 



Random UCF simulation 



This subsection describes a method to generate UCF 
traces 8G{x) from a given correlation function F(8x), 
useful in error analyses such as Fig. [3] 

First, F{8x) is computed over a very wide range 
— L ■ ■ ■ L. Next, the discrete fourier transform 
FFT n {F(Sx)} is computed, which is the power spec- 
trum of UCF (by the Wiener-Khinchin theorem) . Fourier 
amplitudes of UCF, G n = (a n + ib n ) y/FFT n {F(8x)}, 
are then generated with random Gaussian numbers a n , 
b n with a 2 j = 6^ = i. Finally, the conductance fluc- 
tuations are given by the inverse transform 8G(x) = 
lFFT(x){G n }. 
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